Transport coefficients of driven granular fluids at moderate volume fractions 
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In a recent publication [Phys. Rev. E 83, 011301 (2011)], VoUmayr-Lee et al. have determined by 
computer simulations the thermal diffusivity and the longitudinal viscosity coefficients of a driven 
granular fluid of hard spheres at intermediate volume fractions. Although they compare their 
simulation results with the predictions of kinetic theory, they use the dilute expressions for the 
driven system and the modified Sonine approximations for the undriven system. The goal here is 
to carry out this comparison by proposing a modified Sonine approximation to the Enskog equation 
for driven granular fluids that leads to a better quantitative agreement. 



PACS numbers: 05.20.Dd, 45.70.Mg, 51.10.+y, 05.60.-k 

In a recent paper, Vollmayr-Lee, Aspelmeier, and Zip- 
pelius [l| have determined the time-delayed correlation 
functions of a homogeneous granular fluid of hard spheres 
at intermediate volume fractions. As usually done in 
computer simulations, the fluid is driven by means of a 
stochastic external force such that the particles are ran- 
domly kicked with a given frequency ,2J. In the steady 
state, where the external energy input due to driving 
compensates for collisional cooling, they measured the 
dynamic structure factor S{q,uj) (g is the wave number 
and w is the angular frequency) for several values of the 
coefhcient of restitution e and volume fractions rj. The 
corresponding best fit of the simulation results of S{q, lo) 
allowed them to identify (for the smallest values of q) the 
thermal diffusivity Dt and the longitudinal viscosity vi 
coefficients. These transport coefficients are defined as 
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where n is the number density, m is the mass of a parti- 
cle, d is the dimensionality of the system, k, is the thermal 
conductivity, Tyshcar is the shear viscosity and C, is the 
bulk viscosity. The authors compared their simulation 
results for Dt and vi with theoretical predictions of ki- 
netic theory (see Table II of Ref . [1] ) and concluded that 
the latter agree well with computer simulations. In this 
paper, Vollmayr-Lee et al. [11 used three different analyt- 
ical results, two of them 0, |1] obtained from the revised 
Enskog theory (RET) and one of them obtained from 
a (simple) kinetic model of the RET [5|. However, the 
kinetic theory expressions for the transport coefficients 
k(?7) and ?yshcar(?/) employed for making the comparison 
(namely, those derived from the true RETjJ, Hand not 
from a simple model of the latter equation Q) corre- 
spond actually to their low-density forms and so they do 



not incorporate the density corrections accounted for by 
the Enskog kinetic theory [y]. The question arises then 
as to whether, and if so to what extent, the conclusions 
drawn from Ref. [l| may be altered when the Enskog ex- 
pressions for Dt and vi are considered. This question is 
addressed here and additionally, a theory based on the 
Enskog equation is proposed for driven systems that leads 
to a better qualitative agreement with simulation data. 
We consider a granular fluid composed of smooth in- 
elastic disks {d = 2) or spheres {d — 3) of mass m and 
diameter a. Collisions are characterized by a (constant) 
coefficient of normal restitution < e < 1, with e = 1 
in the elastic limit. In the hydrodynamic regime, the 
Navier-Stokes constitutive equations of the pressure ten- 
sor Pij and the heat flux q are given by 



Pij = p5ij -Vshcar [ ^]Ui + ViUj - -%V • U ) -(^V-U, 



q = -kVT - ^iVn, 



(3) 
(4) 



* Elec tronic address: |vicenteg@unex.es] 

URL: |http : //mvl . unex ■ es/e«eb/f isteor/vicente/ 1 



where p is the hydrostatic pressure, u is the flow ve- 
locity, and T is the granular temperature. Moreover, 
apart from the shear and bulk viscosities and the ther- 
mal conductivity coefficient, the coefficient /i is a new 
transport coefficient not present in the elastic case. Ex- 
plicit expressions for the above transport coefficients can 
be obtained from kinetic theory. In particular, the RET 
provides a reliable description of transport properties of 
a granular fluid at moderate densities and across a wide 
range of degrees of dissipation. The RET has been solved 
by applying the well-known Chapman- Enskog method [7| 
conveniently adapted to dissipative dynamics. First at- 
tempts |8|-|10| to determine the Enskog transport coef- 
ficients were carried out for nearly elastic particles so 
that the results only hold in principle in the quasielastic 
limit (e < 1). Subsequent works [ll|, [l2| based on the 
Chapman-Enskog method do not impose any constraints 
on the level of dissipation and take into account the (com- 
plete) nonlinear dependence of the transport coefficients 
on the coefficient of restitution e. On the other hand, 
as for elastic collisions [3], the transport coefficients are 
given in terms of the solutions of a set of coupled linear 



integral equations that are solved by means of a polyno- 
mial Sonine expansion. For simplicity, usually only the 
lowest Sonine polynomial (first Sonine approximation) is 
retained [Il|, [I4I and the results obtained from this sim- 
ple approach agrees well with Monte Carlo simulations, 
except at high dissipation for the heat flux transport coef- 
ficients [13]. Motivated by this disagreement, a modified 
version of the first Sonine approximation was proposed by 
Garzo et al. [^. This modified Sonine approach signifi- 
cantly improves the e-dependence of k and /x and corrects 
the discrepancies between simulation and theory for a di- 
lute gas. The theory proposed in Ref. [4] (in its dilute 
version) was employed in the comparison carried out by 
VoUmayr-Lee et al. [i| in Table II of their paper. 

An important point in all the above theories [J, [ill ll3 
is that they have been derived from an expansion around 
the free cooling state, namely, a reference state that de- 
pends on time through its dependence on granular tem- 
perature. On the other hand, as said before, many of the 
computer simulations (as the one carried out in Ref. [l|) 
are performed by using stochastic driving mechanisms 
to maintain the system in a stationary state. Given that 
this external energy input does not playa neutral role on 
the dynamics properties of the fluid [Ijl, the transport 
coefficients in the undriven (unforced) case [J, [ll|, [iJI 
differ from those obtained in the presence of an exter- 
nal thermostat. In particular, the coefficients jyshoan ^, 
K, and /i were determined in the Appendix B of Ref. [3] 
for a dense fluid of hard spheres (d = 3) heated by a 
stochastic thermostat. However, as in Refs. [Ill, [i2|j the 
results were obtained by considering the (standard) first 
Sonine approximation. To the best of my knowledge, the 
corresponding modified version to the driven case has not 
been derived so far. One of the goals of the present paper 
is to implement this modified version of the first Sonine 
approximation when the fluid is driven by means of the 
stochastic thermostat. Given that the derivation of this 
theory follows similar mathematical steps as those made 
in the undriven case [J], here only the final expressions 
for the transport coefficients jyshoan Cj ^-^d k, will be dis- 
played. 



Since the modified Sonine approximations applied to 
the RET for the undriven 4] and driven cases are per- 
haps the two most accurate approaches for the Enskog 
transport coefficients, only both theories are used in this 
paper. Moreover, for the sake of comparison, the same 
units as those taken in Ref. [i| {m = a = T = 1) are 
chosen. In a compact form, the coefficients 77shcar, Ci ^-i^d 
K can be written as 
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In the above equations, rj — [■k'^^'^/2'^ ^dr(d/2)]n is the 
volume fraction and xM is the pair correlation function 



at contact. In addition, h^ — b^ = and 



c(e) 
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in the case of the driven (stochastic thermostat) RET. In 
the case of the undriven RET f'^, one has 6,, = i, &«, — 2, 



C(^-'?) = Z^2^''^-VX(1 + .) 
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FIG. 1: (color online) Plot of the thermal diffusivity Dt as 
a function of the volume fraction rj for two different values 
of the coefficient of restitution e: e = 0.9 (a) and e = 0.8 
(b). The solid and dashed lines are the theoretical predictions 
derived for the driven and undriven gas, respectively. Symbols 
are simulation results obtained by VoUmayr-Lee et al. [1| at 
e = 0.9 (circles) and at e = 0.8 (triangles). 




FIG. 2: (color online) Plot of the longitudinal viscosity vi as 
a function of the volume fraction r; for two different values 
of the coefficient of restitution e: e = 0.9 (a) and e = 0.8 
(b). The solid and dashed lines are the theoretical predictions 
derived for the driven and undriven gas, respectively. Symbols 
are simulation results obtained by VoUmayr-Lee et al. [3| at 
e = 0.9 (circles) and at e = 0.8 (triangles). The dotted lines 
are the theoretical results (blue for a = 0.9 and orange for a = 
0.8) obtained from the driven theory but when one considers 
the low density value of the shear viscosity coefficient. 



and 



d + 2, o, / 3 
4d ^ ^-^ V 32 



(12) 



For hard spheres (d = 3), as in Ref. [IJ, the Carnahan- 
Starhng approximation 15] for x is considered: 



1 



X = 



(1 - r^r 



(13) 



In the three-dimensional case, according to Eqs. ([T]) 
and (O, the thermal diffusivity Dt and the longitu- 



dinal viscosity vi are defined as Dt = (7r/9?7)K and 

l^l = (7r/677)[(4/3)7?,hear + C]- 

Table |I] reports the comparison between molecular dy- 
namics (MD) results and theoretical predictions given by 
the RET in the undriven [J] and driven cases for hard 
spheres. Given that the fits of the simulation results re- 
quire g-dependent transport coefficients, here only the fit 
result for the smallest value of the wave number q is con- 
sidered. An inspection of the results displayed in the ta- 
ble shows that both kinetic theory results (undriven and 
driven RET) reproduce the trends of simulation data, 
except perhaps in the case of the longitudinal viscosity 
1^1 at the largest volume fraction considered (ry = 0.2). 
At a more quantitative level, it is apparent that the an- 
alytical results derived in the presence of the stochastic 
thermostat agree better with MD than those obtained 
in the unforced case. This is the expected result since 
the simulations were carried out under the action of this 
driving mechanism. It is especially significant the good 
agreement found for the thermal diffusivity coefficient 
Dt, where the discrepancies between theory and sim- 
ulation are less than 10%, even for strong dissipation 
{a — 0.8). It must be noted that the results derived for 
Dt in the undriven case from the kinetic model of the 
RET agree better with simulation data than those 
obtained from the true Enskog equation for an undriven 
granular gas|4|. This is essentially because the impact 
of dissipation on Dt predicted by the kinetic model is 
less significant than the one obtained from the RET for 
undriven systems, which agrees qualitatively with MD 
simulations of driven fluids [I| . 

In order to illustrate the complete 77-dependence of the 
above transport coefficients, Figs. [T] and [H show Dt and 
Vi versus 77, respectively, for a = 0.9 and 0.8. As said 
before, we observe that in general the influence of dissi- 
pation on the thermal diffusivity Dt in the driven case 
is more important than in the undriven case, which is 
consistent with MD data. The influence of the ther- 
mostat on the longitudinal viscosity vi is less signifi- 
cant than for Dt since both theories agree very well. 
On the other hand, the driven RET is still superior to 
the undriven RET, although the former overestimates 
the simulation data for moderate densities. Surprisingly, 
the comparison agrees better when one uses the dilute 
shear viscosity form ?7shoar(0) instead of the correspond- 
ing Enskog coefficient rishcariv) fo^' the coefficient 1^1(77). 
This is clearly shown in Fig. [2] where the combination 
(7r/677)[(4/3)ryshear(0) -I- Civ)] is also plotted for compar- 
ison. The disagreement between the driven RET and 
MD at moderate densities could be due to the fact that 
the density dependence of the shear viscosity is not well 
captured by the modified Sonine solution to the Enskog 
equation or this can also reflect the limitations of the 
RET as the granular fluid becomes denser. In this latter 
case (strong dissipation and moderate densities) , velocity 
correlations among the particles which are about to col- 
lide (which are absent in the Enskog description) could 
play a significant role in the dynamics of the system j]J| . 



On the other hand, more comparisons between the results 
derived from the Enskog equation and computer simula- 
tions are needed before qualitative conclusions can be 
drawn. 

To summarize, I have revisited the comparison for the 
thermal diffusivity Dt and the longitudinal viscosity vi 
carried out in Ref. \V\ between kinetic theory and MD 
simulations for driven granular fluids. As Eqs. ([T]) and 
([2]) show, Dt and vi are defined in terms of the shear 
and bulk viscosities and the thermal conductivity coef- 
ficient. Although the simulations of Ref. [l[ were per- 
formed at moderate densities, in order to compute Dt 
and vi Vollmayr-Lee et aZ.[l|] used the dilute expressions 
of the shear viscosity and the thermal conductivity coeffi- 
cients instead of their corresponding Enskog forms. The 
present comparison (displayed in Table |T] and in Figs. 
[1] and [2]) shows that in general the resulting transport 
coefficients obtained from the RET in the driven case 
(by using a modified Sonine approximation) agree better 
with computer simulations (which were carried out by 
introducing a stochastic driving mechanism) than those 
derived in the undriven case [Jl- This contrasts with the 
comparison made in Ref. [l| where the results obtained 
for an unforced gas were shown to be superior to the 
other ones. On the other hand, at a more quantitative 
level, while the driven theory compares very well with 
simulation data for Dt in the wide range of densities, 
some discrepancies (see Fig. [2|) appear in the case of vi 
as the gas becomes denser. These discrepancies can be 
mitigated if one considers the dilute form [r] = 0) for 
the shear viscosity r/shcar in the definition ^ oi vi. This 



suggests that perhaps the dependence of ryshoar on the 
density r/ is less important than the one predicted by the 
Enskog equation and consequently, the ry-dependence of 
vi is essentially given by the bulk viscosity ^ [which is 
well captured by its Enskog form ([5|)]. 

Before closing this paper, let me offer some remarks on 
the comparison made here between kinetic theory and 
MD simulations. As in the case of ordinary (elastic) 
fluids, the Enskog kinetic equation takes into account 
spatial correlations through the pair correlation function 
but neglects velocity correlations (molecular chaos hy- 
pothesis). In this context, the Enskog equation can be 
considered as an accurate and practical generalization of 
the Boltzmann equation to finite densities. On the other 
hand, although the molecular chaos assumption can be 
questionable as the density of granular fiuid increases 
[16| . there is still substantial evidence in the literature 
on the validity of the Enskog theory for densities outside 
the Boltzmann limit and values of dissipation beyond the 
quasielastic limit [13] ■ The results presented here give 
support again to the Enskog theory as a reliable basis 
for the description of dynamics across a wide range of 
densities, length scales, and degrees of dissipation. 
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